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An  analytic  theory  for  the  large  (pa.ssive)  scalar  derivative  skewress  for  turbulent 
transport  in  a  mean  scalar  gradient  has  been  devised  using  Lie  Algebraic  methods  to 
solve  the  Hopf  equation.  This  Reynolds  independent  skewness  strongly  contradicts 
Kolmoporov ' s  theory  occurs  in  shear  flows,  which  are  being  studied  numerical.ly , 
emphasizing  the  analogies  between  scalar  and  momentum  transport,  ^ 
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I.  Turbulent  Transport  of  a  Passive  Scalar  and  Momentum 

During  this  final,  4th  year,  no  cost  extension  of  my  AFOSR  grant,  my  research  has 
focused  on  the  fiuctuations  of  a  passive  scalar  in  a  linear  externally  imposed  mean  gradient 
and  is  evolving  towards  the  same  considerations  for  momentum  i.e.,  the  fiuctuations  in 
Reynolds  stress  in  a  boundary  layer.  The  two  problems  are  quantitatively  similar. 

A.  Scalar  Fluctuations 

Within  the  last  several  years,  we  have  refined  our  calculations  for  the  pronounced 
exponential  tails  seen  in  the  one  point  stationary  scalar  P.D.F.  sustained  by  a  linear  gra¬ 
dient  and  acted  upon  by  homogeneous  turbulence.  This  effect  was  predicted  by  Pumir, 
Shraiman  and  myself  and  subsequently  verified  experimentally  by  Warhaft  and  Gollub 
with  their  respective  collaborators. 

Our  current  calculations  have  gone  beyond  the  previous  phenomenology  and  express 
via  path  integrals,  the  probability  that  the  history  averaged  strain  acting  on  a  parcel  of  fluid 
is  small.  Such  parcels  can  then  be  transported  from  further  away  along  the  mean  gradient 
and  thus  give  rise  to  tails  in  the  scalar  distribution,  which  we  show  are  exponential.  The 
history  averaged  strain  rate  is  not  a  new  concept  and  has  been  mentioned  in  various  papers 
of  John  Lumley.  Lacking  was  a  means  for  its  calculation,  and  path  integrals  naturally  do 
the  job.  An  essential  ingredient  to  the  observed  exponential  shape  of  the  PDF’s  is  the 
exponential  divergence  of  near  by  trajectories.  This  property  of  non  zero  strain  rate  -  or 
positive  Lyapanov  exponent  is  virtually  a  definition  of  turbulence  and  is  not  shared  by 
pure  shear  i.e.,  Vx  =  ci{t)  +  b{t)y  even  if  time  dependent.  Majda  has  analyzed  this  case  in 
a  number  of  papers  (e.g.,  Phys.  Fluids  A5  1963,  (1993))  and  we  consider  this  velocity  field 
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to  be  a  different  problem. 


The  PDF  of  scalar  differences,  structure  functions,  or  moments  of  the  scalar  gradient 
are  logically  problems  that  follow  a  study  of  the  scalar  itself.  However  some  puzzling  exper¬ 
iments  on  shear  flows  from  the  late  70’s  make  these  topics  quite  central  to  an  understanding 
of  turbulence.  The  experiments  demonstrated  a  Reynolds’  number  independent  derivative 
skewness  or  roughly  linear  3rd  order  structure  function.  These  are  properties  normally 
associated  with  the  longitudinal  velocity  and  any  Kolmogorov  like  theory  for  the  scalar 
predicts  them  to  vanish  strongly  with  Reynolds’  number.  This  effect  is  not  a  particularity 
of  shear  flow  as  recent  numerical  and  laboratory  experiments  at  Cornell  on  homogeneous 
turbulence  show.  However  the  effect  does  disappear  for  a  white  noise  velocity  field. 

We  have  found  path  integrals  ill  suited  to  handling  a  velocity  field  which  scales  and 
has  a  physically  sensible  correlation  time.  Instead  we  employ  the  Hopf  equation  for  the 
three  and  four  point  correlation  functions.  With  a  plausible  representation  for  non  white 
velocity  fields,  the  Hopf  equation  can  be  solved  by  Lie  Algebraic  methods,  so  the  problem 
is  integrable  in  a  sense. 

B.  Shear  Flows 

There  is  no  need  to  emphasize  the  relevance  of  turbulent  boundary  layers  to  all  manner 
of  engineering  and  defense  related  problems.  The  fluctuations  well  into  the  log  layer 
resemble  very  much  the  slightly  simpler  problem  of  homogeneous  shear  as  experimental 
investigations  of  Corrsin  and  simulations  of  Moin  and  others  have  shown.  There  is  no  wall 
in  homogeneous  shear,  but  merely  a  mean  velocity  =  sy  and  an  initial  turbulence  scale 
L  which  grows  down  stream  in  the  experiments  or  in  time  in  the  simulations. 
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Alain  Pumir  during  a  recent  visit  to  Cornell  funded  by  the  A.F.O.S.R.  has  simulated 
a  homogeneous  shear  under  stationary  conditions  with  the  computational  box  setting  and 
fixing  the  integral  scale  of  the  turbulence.  This  had  not  been  done  before.  Naturally  one 
asks  if  the  small  scales,  to  the  extent  one  can  given  the  scale  range,  look  Kolmogorov  like. 
The  velocity  spectra  are  reasonable,  but  the  skewness  Sk  =<  {dy  >  /  <  [dy  v xY  >  = 
is  order  1  and  Reynolds  independent.  A  similar  result  was  found  by  Wallace  experimentally 
in  a  boundary  layer.  The  Vx  involved  is  the  fluctuating  component,  but  with  or  without 
the  mean  velocity  included,  Kolmogorov  theory  predicts  Sk  ^  sj  <  {dyVxY  >^~  Res . 

The  nature  of  this  violation  of  Kolmogorov  theory  is  precisely  analogous  to  the  scalar 
derivative  skewness  and  we  believe  has  a  similar  origin  phenomenologically.  Mixing  of 
scalar  or  Vx  occurs  efficiently  in  blobs  -  making  the  total  Vx  roughly  y  independent  and 
expelling  the  imposed  mean  gradient  into  vortex  sheets  between  the  blobs.  The  blobs  are 
of  course  transient.  Pumir  has  also  found  numerically  that  the  PDF  of  the  Reynolds’  stress 
conditioned  upon  Vx,  Vy  is  very  similar  to  the  Nusseult  number  conditioned  upon  Vy  and 
the  scalar. 

Therefore  the  quintessentially  nonlinear  problem  of  boundary  layer  turbulence  has  sev¬ 
eral  features  in  common  with  scalar  transport.  Further  quantitative  comparisons  should  ^ 
be  possible  when  the  theory  for  the  scalar  has  advanced  further.  Given  such  close  corre¬ 
spondence  between  the  transport  of  scalar  and  momentum  down  a  mean  gradient  there  is 
hope  for  understanding  how  this  similarity  arises  at  the  level  of  the  governing  equations. 
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